Limitations in cooling electrons by normal metal - superconductor tunnel junctions 
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We demonstrate both theoretically and experimentally two limiting factors in cooling electrons 
using biased tunnel junctions to extract heat from a normal metal into a superconductor. Firstly, 
when the injection rate of electrons exceeds the internal relaxation rate in the metal to be cooled, 
the electrons do no more obey the Fermi-Dirac distribution, and the concept of temperature cannot 
be applied as such. Secondly, at low bath temperatures, states within the gap induce anomalous 
heating and yield a theoretical limit of the achievable minimum temperature. 
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Refrigerators are generally characterised by their cool- 
ing power, coefficient of performance, and operating tem- 
perature under various working conditions. To assign a 
temperature to a system, one needs to assume that the 
energy relaxation within the system is faster than any 
rate associated with the heat flux between the system in 
concern and its surroundings. If this condition fails, the 
energy distribution of the particles of which the system is 
formed is non-thermal, and applying the concept of tem- 
perature is strictly speaking inappropriate. Such a limit 
can be achieved in submicron-size coolers at low temper- 
atures. The structure we study is a symmetric configu- 
ration of a NIS refrigerator formed by a series con- 
nection of two Superconductor (S) - Insulator (I) - Nor- 
mal metal (N) tunnel junctions sharing the N island to 
be cooled in between them (SINIS) We demonstrate 
two striking phenomena occuring in these electron micro- 
coolers at low temperatures: evidence of non-thermal en- 
ergy distribution of the cooled electrons and re-entrant 
behaviour with anomalous heating at low bias voltages. 
The cooler performance is typically limited by coupling of 
the electrons to the underlying lattice (electron-phonon, 
e-p coupling). This has, however, strong dependence on 
temperature T: relaxation rate T~_}p slows down on low- 
ering T typically as T~_}p esc Q . Consequently, at low 
enough temperature, characteristically around 100 niK, 
the behaviour of the cooler can be described as if the 
lattice would not exist at all. The interplay of the rates 
for e-p, electron-electron (e-e) and the injection through 
the junctions determines the distribution in the normal 
metal. If e-p or e-e rates are fast, the system assumes 
Fermi-Dirac energy [E] distribution fo{E,Te): 




1 + exp(£/fcBre 



(1) 



In the limit of very strong e-p relaxation, Te equals the 
temperature of the lattice. We call this equilibrium. On 
the other hand, if e-p relaxation is slow, and e-e relax- 



FIG. 1: Scanning electron micrograph of a typical cooler sam- 
ple in (a), and cooling data in (b), where voltage Vp across the 
probe junctions in a constant current bias (28 pA) is shown 
against voltage Vc across the two injection junctions. Cryo- 
stat temperature, corresponding to the electron temperature 
on the N island at Vc = is indicated on the right vertical 
axis. Yet below 100 mK this correspondence is uncertain, be- 
cause of the lack of calibration and several competing effects 
to be discussed in the text. 



ation is fast, Tg is, in general, different from the lattice 
temperature (quasi-equilibrium) . Finally, the slow e-e re- 
laxation as compared to the injection rates implies that 
the electrons assume a non-equilibrium energy distribu- 
tion f{E), which cannot be written as Eq. CQ) U,|5|, and 
it is not possible to assign a true temperature to them. 

Our discussion is motivated by a puzzling experimen- 
tal observation in many coolers at temperatures below or 
around 100 mK Figure 1 shows in (a) a typical SI- 
NIS cooler. The device has been fabricated by standard 
electron beam lithography. The central part forms the 
N island of copper (purity nominally 99.9999 %). The 
injecting tunnel junctions, with normal-state resistances 
R-j (100 fl - 2 k51), contact the island symmetrically into 
the two superconducting aluminium reservoirs at the two 
ends of the N island. The overlapping extra copper shad- 
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ows outside the SINIS structure provide better thermal- 
isation of the Al reservoirs. Two additional NIS probe 
junctions, with normal-state resistances ^ 1 kfi in the 
centre are used to measure the temperature of the N is- 
land or to probe the distribution through the differential 
conductance of the nominally symmetric series connec- 
tion. i?T — 200 n, and thickness of the copper film ~ 30 
nm are the essential parameters of the sample (SI) whose 
data are presented in this paper. The data in Fig. 1 (b) 
were taken by applying a constant current Iq through the 
probe, which, due to the thermal rounding of the current- 
voltage characteristics, provides, by detecting voltage Vp 
across it, a measure of the electron temperature on the N 
island Q ■ We calibrated this dependence by varying the 
bath temperature of the cryostat. Typically one applies 
Iq such that the voltage remains within the gap region of 
the superconductor, Vp < 2A/e, whereby excess heating 
or cooling by the probes is avoided. Here A is the energy 
gap of the superconductor and e is the electron charge 
(A/e ~ 0.2 mV). Figure 1 (b) shows Vp against injection 
voltage Vc; the several curves represent different bath 
temperatures. The approximate temperature calibration 
is given on the right vertical axis. At all but the lowest 
bath temperature of about 50 mK, the curves show the 
expected refrigeration behaviour symmetrically around 
Vn = 0, with optimum cooling at about Vc — ±2 A/e 
[21 • The lowest-temperature curve demonstrates, how- 
ever, a feature which appears as heating at low values 
of Vc and re-entrant cooling again close to 2A/e. This 
behaviour is common with many similar samples, and 
it appears only below 200 mK. Moreover, the measured 
conductance curves of the probe junctions (Fig. 4) indi- 
cate that at low temperatures the actual temperature we 
assign in a measurement depends on the choice of Iq- 
To understand the observed behaviour, we consider the 
properties of the kinetic equation describing the domi- 
nant processes in our system. We assume that the two 
reservoirs are identical and the quasiparticles have a ther- 
mal distribution of Eq. with = Ts on them. In 
steady state f{E) on the N island is then determined by 



[niEp)[fo{En,Ts)-f{E)] 



+n{Ei^)[fo{Ei^,Ts) - f{E)]] = Icon[f;E]. (2) 

Here /coii[/; E] is the collision integral discussed below, 6 
is the level spacing on the island, £'l,r = E zL eVc/2 are 
energies on the left (L) and right (R) reservoirs, and 



n{E) = |Re( 



E + iV 



y^{E + ir)^ - 



)l 



(3) 



is the broadened BCS density of states (DOS) of the su- 
perconductor, r smears the DOS singularity at E = ±A 
and allows for states within the gap, e.g., due to inelastic 
electron scattering in the superconductor or by inverse 
proximity effect due to nearby normal metals. A more 
phenomenological choice is to add a non-zero constant 
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FIG. 2: Energy distribution j{E) against i5/A in (a), (c) and 
(e), and differential conductance G against probe voltage up 
in (b) , (d) and (f ) at three different values of e-e collision pa- 
rameter Jf, which present very slow, intermediate and very 
fast e-e relaxation from top to bottom, respectively. Param- 
eter values r\ =\ - 10~* and Ts/Tq = 0.1 have been assumed. 



F/A to the ideal singular DOS (F = 0). This choice 
would not essentially affect our conclusions. 

It is illustrative to investigate first the case where 
relaxation (both e-e and e-p) is very weak, i.e., when 
Icoi\[f]E] = in Eq. Then we obtain an explicit 

expression for f{E) as 



f{E) 



n{Ep)fo {Ep, Ts) + n{Ei,)fo [E^^ , Tg) 
n{Ep)+n{Ei^) 



(4) 



The solution of Eq. Q is plotted in Fig. 2 (a) for 
five different values of vq = eVc/A, assuming 77 = 
F/A = 1 • 10-4 [SI and Ts = O.lTc- The BCS rela- 
tion A ~ 1.764fcBTc has been assumed for the critical 
temperature Tq. This solution exhibits some nontrivial 
features. At low values of the bias voltage up to z;c ~ 1, 
the distribution first broadens whereafter it starts to get 
narrower again, and at uc = 2.0, it becomes very narrow, 
effectively demonstrating strong cooling. In each case, 
except at t;c = (equilibrium), the distribution is not 
thermal. At v q > 2.0 the distribution would become even 
more unusual [llj, but we do not consider this regime in 
detail for reasons to be explained below. Yet, at those 
voltages the distribution effectively broadens again, sug- 
gesting already the non-monotonic behaviour, heating - 
cooling - heating upon increasing t;c ■ The distribution is 
tested by measuring the differential conductance dl / dVp 
of the probe junctions (Fig. 1 (a)). The relation between 
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dl/dVp, Vp, and f{E) is given by 



(5) 



Here dl/dVp has been scaled by the normal-state con- 
ductance of the series connection of the probe junctions. 
G = has been plotted in Fig. 2 (b) for the colli- 
sionless distributions of Fig. 2 (a) against wp = eVp/A. 
These curves exhibit some similarities to those of biased 
diffusive metal wires 0, especially at low values of vq, 
where quasi-constant DOS at low energies {n{E) ~ F/A) 
mimics resistive normal-metal wire. Another important 
limit is quasi-equilibrium with electrons perfectly decou- 
pled from the phonon bath. The latter condition is well 
justified at low Tg by a standard estimate of the min- 
imum temperature Tj^in determined by the balance be- 
tween cooling power and e-p coupling only 12]. In quasi- 
equilibrium fast e-e relaxation forces the distribution into 
a thermal one (Eq. (Q). Any temperature satisfies this 
and we need to determine T by setting the net heat flux 
from the island, P{T,Ts), equal to zero. In our case we 
can write 



P{T,Ts) 



n{E)[fo{En,T) - fo{E ,Ts)]Ej,dE = 0, 

(6) 

where factor 2 on the right hand side takes into account 
the flux through the two identical junctions L and R. 
The solutions of Eq. © allow us to plot the distribu- 
tions and differential conductance at different values of 
vc as in Fig. 2 (e) and (f). The re-entrant heating- 
cooling-heating behaviour survives but less pronounced 
than in the coUisionless case. It transforms into more 
conventional cooling-heating behaviour above Ts — T*, 
given by 



(A/fcBr*)3/2 exp(-A/fcBr*) ~ ?7/\/27r. 



(7) 



Thus T* depends approximately logarithmically on 77 and 
assumes a value T* /Tq ~ 0.125 when 77 = 1 • 10"*^. 
Equation Q is obtained by equating at low bias the 
quadratic in vc heating and cooling terms arising from 
the states inside and outside the gap, respectively. For 
an intermediate strength of e-e interaction, the distribu- 
tion functions were obtained by Eq. ||2Jl with the e-e 
collision integral given in Ref. '4 . In the case of a dif- 
fusive normal-metal island whose transverse dimensions 
are smaller than the coherence length ^0 = ■\/fi.£'/A, one 
finds /coll — ka/A/coH where 



/coll ~ 



duode 

/(iJ-^A)/((6 + ^^)A)(l 



[/(i5)/(eA)(l - f{E - ^A))(l - /((e -I- a;)A)) 
f{E)){l~f{em 



(8) 
and 



is the dimensionless collision integral, n — — 
D and L are the diffusion coefficient and the length of 
the island. Dividing Eq. Q by (5/ (e^-Rx), we obtain a di- 
mensionless equation where the strength of e-e scattering 
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FIG. 3: Effective electron temperature in tlie cooler against 
the injection voltage. The dashed lines correspond to K = 
(extreme non-equilibrium), K — 0.1 and K = 1.0, from top 
to bottom, and the solid line to K = cx (quasi-equilibrium), 
all with Ts/Tc = 0.1 and r? = 1 • 10"*. The dash-dotted line 
is the result for quasi-equilibrium but for r) — 1 ■ 10~^. The 
horizontal dashed lines indicate the crossover temperature T* 
of re-entrant behaviour in quasi-equilibrium for 77 = 1 • 10~* 
and 1 • 10"'^. The inset shows the dependence of the ultimate 
minimum temperature of the cooler against 77. 



is governed hy K = 2\/2f^y^ IS]. Here i?K = h/e'^ 

and i?T = hD/L^ are the resistance quantum and the 
Thouless energy of the island, respectively. 

Although temperature is not a valid concept in non- 
equilibrium, we can define an effective temperature Toff- 
A natural choice is to require that T^g satisfies the stan- 
dard relation between the temperature and the ther- 
mal energy density of electrons in Sommerfeld expansion, 
which yields 




[f{E) - 1 + 0{E)]EdE. (9) 



Here d{E) is the Heaviside step function. This TcS co- 
incides with the true temperature in (quasi-)equilibrium, 
and it is not affected by the strength of e-e relaxation as 
such. Yet in a biased SINIS, Toff depends on the strength 
of e-e relaxation, because of the heat exchange with reser- 
voirs with non-constant DOS. Figure 3 shows Tes as a 
function of vc at Ts/Tc = 0.1 and = 1 • 10"* for differ- 
ent rates K. In the coUisionless limit, the rise of Tcff is 
largest (almost threefold) and the maximum is reached at 
Vc — 0.8. On increasing the collision rate, the maximum 
Tcff gets lower and it is reached at a lower value of vc , and 
ultimately in quasi-equilibrium (K = 00), the maximum 
is reached at uc ~ 0.4 and its value is about 0.12Tc. The 
minimum temperature in quasi-equilibrium at — 2 
is given approximately by Tmin/Tc ~ 2.57;^/'^ (inset of 
Fig. 3). Although the experimental working tempera- 
tures well below 100 mK imply that 77 < 0.01, we cannot 
make a direct comparison to this theoretical result in the 
present device. For this we would need a thermometer 
calibration down to the lowest temperatures, the electron 
system should be forced to quasi-equilibrium (magnetic 
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FIG. 4: Measured differential conductance of SI at Ts = 340 
mK in (a) , and at Ts — 100 mK in (b) . Zoom-up of a few data- 
sets of (b) are shown in (c) , and the corresponding theoretical 
lines with K = 0.05 in (d). 



impurities would help), and the S reservoirs should be 
well thermalised. 

To assess the degree of non-equilibrium, we measured 
dl/dVp at various values of Vc, which allows for a semi- 
quantitative comparison between the experiment and 
theory (see Eq. lO). In the data, especially below about 
200 mK, we concentrate on the values of both Vc and Vp 
below the gap (< 2A/e) because of the excessive heating 
of the reservoirs when current increases abruptly at the 
gap edge. The data at Ts = 340 mK in Fig. 4 (a) exhibit 
cooling behaviour in quasi-equilibrium: the conductance 
curves measured get first sharper monotonically on in- 



creasing Vc from to 0.35 mV, whereas data at Vc = 0.45 
mV are more smeared, i.e. hotter than any other curve. 
All data sets conform in shape without crossing, demon- 
strating near-to-quasi-equilibrium behaviour. The data 
in Fig. 4 (b), taken at the base temperature of the cryo- 
stat of about 50 mK (best fit to dl/dVp would yield 
Ts ~ 100 mK), show, on the contrary, that the energy 
distribution in this case deviates from the thermal one 
when applying bias Vc- Data at dl/dVp < 0.03 first 
indicate that the low-bias conductance becomes larger 
when increasing Vc from up to 0.1 mV (unlike at 340 
mK), where-after the curves start to get sharper, but 
they heavily cross each other in this regime. At 0.45 mV 
the data present significant heating again. The mere fact 
that the data-sets corresponding to different values of Vc 
criss-cross in the regime below 0.4 mV is a demonstration 
of non-equilibrium. A few data sets from (b) are mag- 
nified in (c), and the corresponding theoretical results, 
assuming K — 0.05, have been shown in (d). The re- 
semblance is obvious, although the theoretical lines show 
abrupt rise from dl/dVp = at Vp = Vc due to the in- 
fluence of the gap edge. This feature is smeared in exper- 
iment most likely because of noise and finite excitation 
level (25 /iV p-p) in the measurement. 

In summary, we have shown that slow electron-electron 
relaxation restricts the use of the concept of temperature 
in electron coolers at low temperatures, and that the non- 
zero DOS within the gap of the superconducting reser- 
voirs gives rise to anomalous heating and determines the 
ultimate minimum temperature that can be achieved. 
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